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ABSTRACT 



Context. Small-scale magnetic fields in the solar photosphere can be identified in high-resolution magnetograms or in the G-band 
as magnetic bright points (MBPs). Rapid motions of these fields can cause magneto-hydrodynamical waves and can also lead to 
nanoflares by magnetic field braiding and twisting. The MBP velocity distribution is a crucial parameter for estimating the amplitudes 
of those waves and the amount of energy they can contribute to coronal heating. 

Aims. The velocity and lifetime distributions of MBPs are derived from solar G-band images of a quiet sun region acquired by the 
Hinode/SOT instrument with different temporal and spatial sampling rates. 

Methods. We developed an automatic segmentation, identification and tracking algorithm to analyse G-Band image sequences to 
obtain the lifetime and velocity distributions of MBPs. The influence of temporal/spatial sampling rates on these distributions is 
studied and used to correct the obtained lifetimes and velocity distributions for these digitalisation effects. 

Results. After the correction of algorithm effects, we obtained a mean MBP lifetime of (2.50 ± 0.05) min and mean MBP velocities, 
depending on smoothing processes, in the range of (1 - 2) kms~'. Corrected for temporal sampling effects, we obtained for the 
effective velocity distribution a Rayleigh function with a coefficient of (1.62 ± 0.05) km s -1 . The x- and y- components of the velocity 
distributions are Gaussians. The lifetime distribution can be fitted by an exponential function. 

Key words. Sun: photosphere, Sun: magnetic fields, Techniques: image processing 



1. Introduction 

Magnetic bright points (MBPs) are small-scale magnetic fea- 
tures, visible as bright points in the solar photosphere. Located 
in intergranular lanes, they are called "bright points" (BP), "net- 
work bright points" (NBP), "G-band bright points" when ob- 
served with a G-band filter, or filigree in a network region. 
MBPs are called "flowers" i f they appear group ed in a roundly 
shaped formation (see e.g. IBerger et al.l I2004). The reported 
diameters of MBPs range from 100 up to 300 km (se e e.g. 
iMuller & Keii 1 19831 IWiehr et al.l 12004 lUtz et all l2009al) . and 
the corresponding magnetic field strength reac hes values of sev- 
eral kG (as revealed by inversion techniques. iBeck et al.l 12007, 
IViticchie et al.l |2009|). MBPs appe ar at the merging points of 
granules (IMuller & Roudien [l992) and display a complex evo- 
lution. Pushed by granules, they are able to form groups, but 
merging and splitting of single features frequently occur as well 
dBerger & Titldll996b . 

It is well known that the MBP brig htness in the G-band is in 
close relation with magnetic fields (lKellerlll992llYi & Engvoldl 
1993llBerger& Title 2001 , B erger et al.l2004lBharti et al.l2 006. 
Viticchie et al.l l2009). This part of the electromagnetic spectrum 
comprises CH molecule lines centered around A = 430 nm. The 
increased brightness of magnetic acti ve features in the G-band is 
due to a decreased opacity (see e.g. lRutten|[T999l ISteiner et al.l 
1200 U ISchtissler et atll2003h . G-band images yield a better res- 
olution than magnetograms, and thus provide good means to 
detect small features such as MBPs. However, one must be 



cautious, because also non-magnetic b rightenings can be f ound 
(e. g. on top of granular fragments, see Berger & Title 2001). 

MBPs are described either by dynamical or statical MHD 
theories. Earl y works on theoret i cal aspects were done e.g. b\ 



ISpruitl d 19761)". lOsherovich et al.1 (1 19831) . iDeinzer et al.l dl984al 



and Ipeinzer et a" | (Q984b). Simulations of dynamical pro - 
cesses in the photo sphere were done by ISteiner et al.1 d 1998b . 
Shelv aget al.l (120061) and bv lHasan et al.l d2005l) . who focused on 
dynamical chromospheric phenomena triggered by photospheric 
flux tubes. 

MBPs are in many ways of interest. First of all it is sup- 
posed that rapid MBP footpoint motion can excite MHD waves. 
These waves could contribute in a significant way to the heating 
of the solar corona (see e.g. IChoudhuri et ailll993l IMuller et al. 
11994 . In addition to the wave heating processes also na noflares 
can be trig gered by these flux tube motions (see e.g. iParketj 
fT98l [19881 . Another research field in solar physics deals with 
the generation of the magnetic fields. Are small-scale fields like 
MBPs connected to the global magnetic field, or are they gen- 
erated by small-scale magneto-convection (see e.g. lVogler et al.l 
2005)? Does the magnetic field interact on small scales with the 
solar granulation pattern and hence influence the solar irradia- 
tion and/or activity? 

In this work we concentrate on the lifetime a nd velocity dis - 
tributions of MBPs like several authors before: IMuller! (1983), 
IMuller et al.1 (1 19941) . IBerger & Title] (119961) . iMostl et al.1 d2006l) . 
In addition to remeasuring these parameters with seeing-free ob- 
servations from the 50 cm Hinode/SOT space-based telescope 
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in the quiet sun, we investigate the dependence of the mea- 
surements on observational aspects like temporal and spatial 
sampling. This is important to understand why different studies 
yielded different results and how these differences could be ex- 
plained. In Sect. 2 we give an overview of the data we used and 
on the stability of the Hinode satellite pointing. In Sect. 3 we de- 
scribe the developed tracking algorithm. In Sect. 4 we present 
our results and the implications of different spatial/temporal 
sampling rates and a way to correct for digitalisation and algo- 
rithm effects. In Sect. 5 we discuss our results in the context of 
other studies, and in Sect. 6 we give a brief summary and our 
conclusions. 

2. Data 

We used two different data sets near the solar disc centre ob- 
tained by the solar optical telescope (SOT; for a description see 
llchimoto et alJl2004ISuematsu et alll2008l) on board the Hinode 
satellite dKosugi et al.ll2.Q07b . SOT has a 50 cm primary mirror 
limiting the spatial resolution by diffraction to about 0.2 arcsec. 
Data set I, taken on March 10, 2007, has a field-of-view (FOV) 
of 55.8 arcsec by 111.6 arcsec with a spatial sampling of 0.108 
arcsec per pixel. The complete time series covers a period of ap- 
proximately 5 hours and 40 minutes, consisting of 645 G-band 
images recorded with a temporal resolution of about 32 seconds. 
Data set II, from February 19, 2007, has a FOV of 27.7 arcsec by 
27.7 arcsec with a spatial sampling of 0.054 arcsec per pixel and 
consists of 756 exposures over a period of roughly 2 hours and 
20 minutes. Data set II has a better temporal resolution of about 
1 1 seconds. Both data sets cover a quiet sun FOV and were fully 
calibrated and reduced by Hinode standard data reduction algo- 
rithms distributed under solar software (SSW). 

We focus now on the Hinode/SOT image stability. As MBP 
velocities are in the range of a few kilometers per second, the 
stability of the satellite pointing and the resulting image stabil- 
ity has to be checked very carefully. This was done by calcu- 
lating the offsets between each image and the succeeding one 
by cross-correlation techniques (e.g. SSW routine cross_corr). 
FigureQ]illustrates the outcome. The calculated offsets between 
two frames were on average lower than one pixel (0.2 + 0.1 pix- 
els). The offsets versus time is plotted in the second row of Fig. 
Q] The se agree with the SOT specified and measured image sta- 
bility (Ichimot o et al.ll2008l) . On the other hand, it is not only 
necessary that the images are coaligned from time step to time 
step, but there should also be no systematic drift of the satellite- 
pointing. In order to check this, we calculated the cumulated dis- 
placements. This showed that the satellite drifted away from the 
original pointing (see top and bottom row of Fig.Q]) with a mean 
speed of about 0.4 km s _1 , which lies in the range of MBP ve- 
locities. Therefore, we had to correct for this drift in the MBP 
data analysis. 

3. MBP tracking algorithm 

For the segmentation and ident i fication of MBPs we used the al- 
gorithm described in lUtz et ail (|2009a). This algorithm was fur- 
ther extended to track identified MBPs in subsequent images, as 
described below. 

After the application of the segmentation and identification 
routines dUtz et alj2 009a). we obtain a set of single analysed im- 
ages. Analysed images means that we have identified the MBP 
features and obtained relevant parameters like size, brightness 
and position. For dynamical parameters like lifetimes and ve- 
locities, it is necessary to analyse not only a single image but a 
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Fig. 1. Hinode/SOT pointing stability for data sets I (left) and 
II (right). Top row: the displacement between an image and the 
succeeding one (Ay versus Ax). The dotted lines indicate sta- 
ble satellite pointing. It can be seen that the actual barycen- 
ter of the displacements is shifted diagonally out of this point 
(drift). Middle row: effective ( -y/Ax 2 + Ay 2 ) displacements ver- 
sus time. Bottom row: cumulated displacements in x- and y- 
direction, which reflects the satellite pointing on the solar sur- 
face. 



series of images. The basic idea behind our tracking algorithm 
is simple. The algorithm takes the position of a detected MBP 
and compares it to the positions of MBPs in the succeeding im- 
age. If two MBPs match within a certain spatial range (which 
was set to 2.8 Mm for both data sets), the points are identified 
as the same MBP at different times. Therefore, a subsequent ap- 
plication of this comparison (for all images and all MBPs in the 
images) leads to complete time-series of MBP features. The al- 
gorithm is able to discern between isolated MBPs and groups of 
MBPs (this is explained in more detail in Sect. 4). The tracking 
algorithm has to consider several aspects as for instance noise 
and finite size of the images as well as the time span of the im- 
age sequence, to work correctly. All features which are identified 
either in a single image or in two images should be regarded as 
noise (or identification artefacts) and are therefore disc arded (for 
possible causes of noise features see lUtz et al.ll2009ah . The spa- 
tial limitation of images can cause the time series of a MBP to be 
interrupted and splitted up. This occurs e.g. when a MBP moves 
across the border of the image. Thus, all MBPs that move too 
close to an image border are dismissed. The beginning and end 
of a time series should be treated the same way. As we do not 
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Fig. 3. Histogram of the obtained MBP lifetimes of data set I 
together with an exponential fit (solid line). The histogram was 
normalised by dividing the counts in each bin by the binsize. The 
histogram is only shown up to the first zero crossing, i.e. there 
are a few longer-living MBPs than shown in the figure. 
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Fig. 2. Top panel: Starting frame of one of the identified MBPs. 
Bottom panels: Evolution (time series) of the detected MBP 
in the rectangle. The white cross marks the derived brightness 
barycenter. The path of the MBP is shown as a white line. 
Positions are given in arcsec from the bottom left corner of the 
full frame. 



know whether a MBP existed already before the first exposure 
was taken, or exists still after the last exposure was taken, all 
MBPs starting with the first image or ending with the last im- 
age are dismissed. All these effects have to be taken into account 
to avoid an underestimation of the obtained lifetimes. Figure [2] 
shows an example of MBP tracking. 
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Fig. 4. Cumulated number of MBPs versus lifetime fitted with 
an exponential function (solid line) for data set I. The values of 
the fit represent the number of MBPs that have a lifetime greater 
than or equal to the corresponding x-axis value (j ^jfdt}. Only 
the statistically robust bins (> 100 data points) were considered 
for the fit. 



4. Results 

4.1. MBP lifetimes 

After identification and tracking of MBPs in the G-band image 
sequences, we derived the lifetime of each MBP. The lifetimes 
were measured by the time difference between the first and last 
detection of each MBP. For a correct interpretation, one should 
distinguish between isolated and grouped MBPs. An MBP that 
does not have any neighbours during its whole lifetime counts 
as an isolated point feature. A grouped MBP is a feature that 
has at least in one instance of its life another MBP in its vicin- 
ity. The vicinity was defined to be approximately 2800 km. In 
this paper, we concentrate on isolated MBPs, as it is quite diffi- 
cult to obtain and interpret the lifetimes of grouped MBPs cor- 
rectly. Grouped MBPs can split up and merge again. Therefore, 
many "substrings" in the group evolution can occur, and it is dif- 
ficult to find a proper definition of the term "lifetime". Various 
kinds of definitions appear possible, e.g. using the average over 
all substrings, the substring with the longest duration, the time 
difference between first and last occurence of a substring. When 
two isolated features meet and form a group, the uniqueness of 
a substring is lost. If they split up again in the further evolu- 
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Fig. 5. Same as Fig. |4]but derived for a reduced sampling rate 
(only each third image was used). The corresponding temporal 
resolution is 96 seconds. The solid line shows a linear fit for the 
statistically more robust data. The dotted line shows a linear fit 
for a second "unresolved" distribution (for details see text). 
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Fig. 6. Cumulated number of MBPs versus time for the two dif- 
ferent data sets and different temporal samplings. Fits belonging 
to data set I are shown by solid lines, data set II by dashed lines. 
The numbers indicate how many images were skipped to artifi- 
cially decrease the temporal sampling. 



tion, it is unclear which one should be followed. In general the 
decision rules will thus strongly influence the derived lifetime 
distribution. Figure[3]shows the measured lifetime histogram for 
data set I with normalised data bins. The number of MBPs with 
a certain lifetime as well as the cumulated (the total (integral) 
number N(t) of all MBPs that have a lifetime of longer than t 
minutes) number of MBPs (see Fig. |4j were fitted with an expo- 
nential function decreasing with lifetime (for an interpretation of 
the fit coefficients see Table[T|i: 

x(t) = a expi-b ■ t). (1) 

The exponentially decreasing number of MBPs can be caused 
by the physics of the process under estimation and/or by the used 
algorithm. If we take into account that the algorithm has a certain 
probability p to find an MBP once, then the probability to find 
the same MBP in the next image again would be p 1 . For n images 
we achieve p" which leads to: 

p{t)=p^\ (2) 

where t = n ■ At is the time from the first MBP detection and At 
is the time difference between two consecutive images. This can 
be easily transformed to: 

p(t)=p-exp^ t Hp))- (3) 

If we suppose that the real number of MBP versus lifetime 
distribution follows an exponential function, A^o exp(-/t ■ t), we 
find for the measured distribution (for a more elaborate deriva- 
tion of the measured distribution see the appendix): 



N{t) = No ■ expi-A p t) ■ p{t) (4) 

= N p ■ sxpi-A p t) ■ ex p(^ m (p)) (5) 

= No ■ p ■ exp + i ln(p)| t j (6) 

= N - p- expi~A tot t). (7) 



In this context A^(f) denotes the measured number of MBPs with 
a lifetime larger than t. No stands for a "fictitious" starting num- 
ber of MBPs. MBPs are formed and decay constantly. If we 
could force all MBPs to be formed in the first exposure and then 
could turn off the forming process, No would denote the number 
of features to be found in the first image. A p gives us the physical 
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Fig. 7. Decaying parameter A, , derived from the fits plotted in 
Fig. [6] as a function of the temporal sampling rate. A first order 
fit was applied, which yields the decaying rate A p corrected for 
detection/algorithm effects (solid line). The derived physical de- 
caying rate value is A p = - 0.40 ± 0.01 min -1 , corresponding 
to a mean MBP lifetime r = 2.50 + 0.05 min. Values obtained 
by data set I are shown as crosses, values derived by data set II 
are plotted by stars. Note that the two data points, derived by dif- 
ferent data sets but with similar sampling rate, closely overlap. 
The 1 -sigma bandwith of the fit is shown by dashed lines, the 
3-sigma bandwith by dash/dotted lines. 

decaying parameter, whereas A, nt is the "total" decaying param- 
eter derived from the measurements with the tracking algorithm. 
A comparision of Eqs. (7) and (6) gives the relationship between 
the measured decaying constant iA tot ) and the temporal sampling 
rate Af as: 

A tot = A p - Inip) (8) 

If we obtain the total decaying constant for data sets with a dif- 
ferent temporal sampling Af, we should be able to recover the 
true physical decaying rate iA p ) as well as the detection proba- 
bility of our algorithm ip). 

To test wheter this relationship really holds, we artificially 
reduced the temporal sampling of our data sets. For each of these 
realisations, we got a distribution similar to Fig.|4] We artificially 
decreased data set I by interleaving one or two images (corre- 
sponds to images every 64 and 96 seconds, respectively) and data 
set II by interleaving one, two and three images (corresponds to 
images every 22, 33 and 44 seconds, respectively). Interestingly 
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Table 1. Explanation of the different fitting parameters used for Eq.Q] 
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distribution type fit parameter representation meaning 



lifetimes 


x{t) 


A(t) 


Number of MBPs with a lifetime of t minutes in a 1 minute interval. 
Could be interpreted as "activity" (according to radio-activity): 
number of MBPs that decay within 1 minute at time t. 




a 
b 


A(0) 
A 


Starting activity (number of MBPs which will decay within the first min) 
Decaying parameter (related to the mean lifetime t = 1 /A) 



cumulated x(f) N(t) Number of MBPs with a lifetime greater than t 

lifetimes a N(Q) Number of all MBPs 

b A Decaying parameter 



Velocity (v x ) histogram of MBPs 
0.5 E 




-6 -4 -2 2 4 

Velocity [km s"'] 



Fig. 8. The top left panel shows the x-component of the velocity of identified single MBPs together with a Gaussian fit for data 
set I. The fit parameters are: [i — - 0.02 km s and cr = 0.89 km s . This distribution was derived by skipping two images 
(giving a temporal resolution of 96 seconds). The y-component is shown below. The fit parameters are: fi = - 0.04 km s and 
cr = 0.85 km s~ . The probability density is displayed on the right side. It can be seen that there is no preferential direction for the 
MBPs movement. 

on thus other realisations is that a tail is forming out of the mea- 
surements cloud of the long-living MBPs of Fig. [4] This can be 
seen e.g. in Fig. [5] which shows the cumulated lifetime for an im- 
age cadence of 96 seconds. This could represent a second (unre- 
solved) distribution, indicated by the dashed line. Since we only 
found about 20 MBPs with lifetimes longer than 15 minutes, we 
did not investigate this tail further. The resulting exponential fits 
to the histograms are shown together with the original ones in 
Fig. [6] It can be seen that the slope of the distributions is very 
similar for skipping images in data set I to skipping 2 images 
in data set II, which relate both to a similar temporal sampling of 
about 30 seconds (32 to 33 seconds; compare also Fig. 0. The 
obtained fit parameters are listed in Table|2] Decreasing the tem- 
poral sampling leads to a flattening of the exponential slope. In 
Fig.|7]we plot the derived decaying constants A tol versus the tem- 
poral resolution, which follows a straight line. From the linear fit 
to these curve we derived A p = -0.40 + O.Olmin -1 , which cor- 



Table 2. Overview of the measured lifetimes of MBPs. 



Temporal cadence 


Fit parameter 


Fit parameter 


At [s] 


No 


A [min- 1 ] 


11 


990 ± 70 


-1.09 ±0.03 


22 


770 ± 40 


-0.76 ± 0.02 


33 


660 ± 30 


-0.62 ± 0.01 


44 


640 ± 90 


-0.56 ± 0.03 


32 


7600 ± 250 


-0.63 ± 0.01 


64 


9300 ± 450 


-0.54 ± 0.01 


96 


8700 ± 750 


-0.45 ±0.01 



responds to a mean MBP lifetime of 2.5 ± 0.05 min. The second 
fit parameter gives us after some transformations the identifica- 
tion probability p of the MBP features in single images. This 
probability has a value of about 90%. 
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Table 3. Overview of the velocity fit coefficients for the Rayleigh 
distributions (shown in Figs. [TOl and [TTTi for both data sets. 



Velocity histogram of MBPs 



Spat, sampling 


Temp, cadence 


Fit parameter 


Mean value 


[arcsec/pixel] 


At [s] 


cr [km s -1 ] 


/i [km s -1 ] 


0.108 


32 


1.12 


1.59 ± 1.20 




64 


0.95 


1.32 ±0.91 




96 


0.88 


1.17 ±0.76 




128 


0.81 


1.08 ±0.67 




160 


0.78 


1.03 ±0.61 




192 


0.72 


0.94 ± 0.56 


0.054 


11 


1.36 


1.77 ± 1.22 




22 


1.14 


1.49 ±0.90 




33 


1.04 


1.35 ±0.84 




44 


1.01 


1.29 ±0.77 




55 


0.98 


1.25 ± 0.73 




66 


0.94 


1.21 ± 0.70 



4.2. Velocities 

The MBP velocity distributions were derived by measuring the 
movement of the MBPs brightness barycenter. The barycenter 
positions were corrected for the satellite pointing drift as dis- 
cussed in Sect. 2. We smoothed the velocities by reducing the 
available image cadence. The velocity distributions can be fitted 
by Gaussian distributions for the x- and y-velocity components 
(see Fig. [8}, respectively. The effective velocity (v = -^/v 
was fitted by a Rayleigh distribution (see Fig. |9j: 



f(v,cr) 



(9) 



In this notation, cr is the standard deviation of the distribution. 

Figure [10] shows the velocity distributions and Rayleigh fits 
for data set I obtained by different temporal sampling rates. For 
detailed information about the fit coefficients derived for differ- 
ent smoothing levels we refer to Table [3] Figure [TT] shows the 
Rayleigh fits of these distributions all together in one plot. It can 
be seen in both figures that the velocity distribution shows a scal- 
ing behaviour with the cr-parameter. As we have only changed 
the temporal sampling rates, we know that there has to be some 
relation between these parameters. We have tried to find this re- 
lation between the cr coefficient of the Rayleigh distribution and 
the corresponding temporal (Af) sampling rate. We applied the 
following empirical relationship: 



cr 1 = cr 1 + 



(10) 



In units of measure D corresponds to a diffusion parameter of an 
ensemble of MBPs. This means that this parameter should de- 
scribe the random movement and spreading of a population of 
MBPs (not to be mixed up with the diffusion of the magnetic 
field of a single MBP). The application of this relation as a fit 
function on our data sets is shown in Fig. Q~2] We obtained the 
fit parameters cr Q = (1.62 + 0.05) km s" 1 and D = (350 ± 20) 
km 2 s-\D agrees wit hin the errors with the diffusion coefficient 
in lBereer et alj d 1998b . who reported a value of 285 km 2 s _1 for 
a quiet Sun FOV. The interesting question is which effective ve- 
locities would be measured, if we could continously observe the 
Sun rather than just observing the positions of the MBPs at cer- 
tain instances of time? Therefore, we can calculate the limit for 
Af — > that will result in <x . <t should therefore be the "true" 
standard deviation for the velocity distribution independent of 
temporal sampling. Interestingly this value is more than twice 




Fig. 9. Effective velocity v = Jv? + v,. 2 of the MBPs of data set 
I is shown together with a Rayleigh fit. The statistic is defined by 
skipping two images between measurements. The fit parameter 
cr is in this case 0.88 km s _1 . The fit coefficients for different 
smoothing levels are summarised in Table [3] 



(2.2) as large as the smallest measured sigma value and still 1 .4 
times as large as the highest measured value (without consider- 
ation of the sigma value for the highest cadence as this value 
seems to be an outlier), see also Table [3] From this extrapolated 
"true" sigma parameter for the Rayleigh distribution we can eas- 
ily derive a new "true" mean velocity by the equation for the 
theoretical mean of this distribution (fi — cr ■ y/(n/2)), giving a 
value of 2.0 km s _1 . Figure [TH shows the Rayleigh distribution 
created with the estimated true cr-parameter. 



5. Summary and Discussion 

We have further extended the detection a nd iden tification algo- 
rithm for MBPs described in lUtz et al.1 (l2009al) . It consists of 
an additional processing step, which enables us to analyse the 
time series of G-band images. We found MBP velocity and life- 
time distributions in accordance with earlier publications. The 
lifetime distribution shows an exponential decrease for longer 
lifetimes. The decaying parameter was found to be 0.4 min -1 
corresponding to a mean lifetime of 2.5 min. We outlined a way 
to estimate the identification probability of the algorithm (about 
90%) and to correct the lifetimes for the effects of imperfect de- 
tection. 

If we compare our lifetime-result with previous studies (see 
Table 3J, we see that lifetimes obtained by automated feature 
tracking algorithms tend to be shorter than those measured by 
visual identification and tracking. The reason for this could be 
that if we identify a feature in a certain exposure visually, we 
would try to find the same feature in the previous and subsequent 
images. Since we presume that in the next or previous image the 
same MBP feature will be found, we tend to identify structures 
as beginning or decaying MBPs (even when they are too weak to 
be identified by an automated identification routine or are pos- 
siblly unrelated features) that we would not identify if we had 
no clue about the temporal evolution. The results also depend 
on how the network and inter-ne twork fields are defined and 
observed. In de Wiin et al.1 d2008l) small-scale fie lds were first 
identifi ed in magnetograms; on the other hand, in de Wii n et ail 
(2005) MBPs were used as proxies for the magnetic field result- 
ing in a difference of the derived lifetimes of a factor of three. 

Our results suggest that the different values for MBP life- 
times are related to different data sets (sampling) and methods 
but may also depend on the amount of magnetic flux in the re- 
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Fig. 10. MBP velocity distribution for different smoothing levels for data set I. The temporal resolution is reduced from the top left 
panel (each image used; At = 32 s) to the bottom right (five images are skipped between successive frames; Af = 192 s). Solid 
curves represent the fits performed by Rayleigh functions. The fit parameters are given in Table[3] 



Tabl e 4. Overview on MBP lifetimes in earlier studies. In 
IdeWijnetaU (2008) magnetograms were used for the identifi- 
cation of the small-scale magnetic fields (instead of MBPs). 



Paper reported MBP lifetime [min] 

de Wiin et al. (2008) 10 

Mostl et al. (2006) 4.4 (±2.4) 

de Wiin et al. (2005) 3.5 

Sanchez Almeida et al. (2004) < 10 

Berger & Title (1996) ~ (6 - 8) 

this study 2.5 



gion of interest, i.e. depending on the local physical properties 
of the solar photosphere. Near large magnetic flux regions such 
as plages or sunspots, MBP lifetime could be larger as a result 
than in quiet Sun regions. This can be explained if one keeps in 
mind that the dissipation of large amounts of flux would in turn 
also take longer. Our re sults, obtained from the analysis of iso- 
lated MBPs agree with Ide Wiin et ail (120051) . In this work, the 



authors excluded network regions (i.e., MBP chains) from the 
analysis through the use of Ca II H maps. The s electi on of iso- 
lated bright features in both the Ide Wiin et alj d2005l) analysis 
and the analysis reported here is the most plausible argument to 
explain the agreement above cited. The second unresolved distri- 
bution at larger lifetimes in Fig. 5 may be related to small patches 
of active network in our datasets. Although we aimed to exclude 
active patches (characterised by grouped MBPs) in our study, 
it could happen that some points were still identified (we note 
that Fig. [5] shows in total only about 20 of these events). The 
derived effective velocity distribution was fitted by a Rayleigh 
distribution. We saw that the velocity distribution depends on 
the temporal resolution. Depending on this sampling effect we 
estimated the Rayleigh fit parameter to be in the range of 0.7 
to 1.4 km s _1 . This corresponds to mean velocities of 0.9 to 1.8 
km s _1 . The true cx-parameter for the Rayleigh velocity distri- 
bution was estimated to be 1.6 km s _1 , which corresponds to a 
mean velocity of about 2.0 km s _1 . 

To estimate the velocity of the MBPs we used the brightness 
barycenter to get their positions and subsequently their veloc- 
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ities. Alternatively we could have measured the movement of 
the brightest pixel of an identified feature. Both methods have 
their advantages and disadvantages. Measuring the movement 
by tracking the brightest pixel of a feature is unique and not de- 
pendant on the shape of a feature. On the other hand this method 
leaves only a few possibilities of movement due to the achiev- 
able spatial sampling (in principle: up/down, left/right, and di- 
agonal; one or two pixels). Therefore, it would be hard to gain 
a "real" velocity distribution. By taking the method of tracking 
the barycenter of brightness, many more possibilities are avail- 
able, as the barycenter coordinates can be determined on sub- 
pixel resolution. However, this method is sensitive to the mor- 
phology of the feature and therefore on the definition of the size 
of the feature (see also lUtz et alj2009bl on different possibilities 
of size definition). 

Our tracking algorithm estimates the path of a MBP by 
comparing the positions of a MBP in two consecutive images. 
Therefore one is obliged to choose a certain comparison range, 
which could be easily estimated by multiplying the maximum 
anticipated velocity (about 4 km s -1 ) by the temporal resolution. 
If the range is set too small, time series of MBPs will be broken 
up, causing the lifetimes and velocities to be underestimated. On 
the other hand, by choosing the range to large, unrelated features 
will be connected to a time series, which leads to an overestima- 
tion of lifetimes and velocities. 

To get a real velocity distribution rather than just a few veloc- 
ity "states'Q we have to smooth our obtained velocity values. The 
smoothing of the distributions has several effects. First of all the 
distribution gets "closer" to the real distribution as there are now 
more possible velocity states. On the other hand the measured 
velocity decreases and the distribution gets narrower. This can 
be explained by the movement of a feature on the solar disc. This 
movement consists of two parts, a deterministic movement and 
a chaotic movement component. The deterministic movement is 
due to large scale flows, e.g., supergranular flows, meridional 
circulation etc.. These can cause MBPs to move over larger spa- 
tial scales and longer time scales with preferred directions and 
velocities than just moving randomly. The chaotic movement is 
due to random processes and corresponds to a "zigzag" path of a 
feature on the solar disc (for implications of t his random walk on 
the diffusive magnetic field transpo rt see e.g. lBerger et alJ ll998. 
IStanislavskv & Weron 2007, 2009). This yields a dependence of 
the length of the path on the temporal resolution of the obser- 
vations. Consequently the shorter the temporal lags between the 
observations, the more accurate is the measurement of the path 
and the longer the path will be (one could observe more zigzag 
steps). So we have learned that by smoothing (artificially de- 
creasing the time cadence), the derived MBP velocity decreases. 
This is not only true for MBP movements but h olds also for 
other motions, as can be seen e.g. lAttie et alJd2009h who investi- 
gated granular flow field maps and found a similar behaviour. On 
the other hand, if we do not smooth, we probably get velocities 
which are too high. This is due to the effect that the temporal to 
spatial sampling has to have a certain ratio. If we have a temporal 
sampling which is too high with respect to the spatial sampling 
(for a given velocity), the MBP would stand still for several ex- 



1 A state is a possible velocity value. As we are working with discre- 
tised data, there are no arbritrary displacement possibilities and there- 
fore also the velocity values are not arbritary. Some velocity values are 
a lot more favourable (no movement, one pixel in a direction,. . . ) than 
others. Calculating displacements by barycenters instead of calculat- 
ing them by the difference of the brightest pixel helps, but as the sizes 
are very small the displacement states mentioned before are still more 
favourable. 



Effective velocity of MBPs 
1.01 




Velocity [km s" 1 ] 

Fig. 11. Comparison of all Rayleigh function fits to the velocity 
distributions plotted in Fig. [10] The broadness of the distribution 
as well as the mean velocity of the MBPs decreases with increas- 
ing smoothing level (from red to blue), i.e. lower time cadence. 



posures. If the feature finally moves, it jumps from one pixel to 
another and the derived velocity will be too high, since in real- 
ity the feature moves more or less continuously all the time. On 
the other hand, if the time resolution is too low (smoothing too 
high), we would loose much of the chaotic movement and de- 
rive a "mean" velocity which is too low. Figure[10]shows this ef- 
fect. We see that in the first cases the distribution is broader with 
a high velocity tail. The last two plots show that the obtained 
measurement points now concentrate more or less at average ve- 
locities (the distribution gets significantly narrower) and only a 
small number of high velocities can be found. 

Another interesting aspect is that not only the shape of 
the distribution but also the mean value of the distribution is 
changed. This is due to the fact that we have a non-symmetric 
function. In the process of smoothing, more of the high val- 
ues are redistributed to lower values than vice versa. This is the 
mathematical explanation. The technical interpretation is that we 
loose a part of the MBPs path by reducing the temporal sam- 
pling. This leads to the shift of the distribution which can be 
seen in Fig.fTTl 

6. Conclusion 

In this study we found MBP lifetimes (corrected mean value 
2.50 + 0.05 min) and velocities (corrected Rayleigh fit param- 
eter 1.62 ± 0.05 km s _1 ) in the range of earlier reported values. 
Additionally we investigated the relationship between image ca- 
dence and obtained lifetimes and velocities. Using the derived 
relations we were able for the first time to correct the measured 
velocity distribution for this influence. The corrected velocity 
distribution shows many more fast moving MBPs (v > 4 kms -1 ) 
than the original measured distribution (12% for the corrected 
distribution; 0.2% for the measured distribution with the highest 
cr; see also Fig.[T3l). The obtained fraction of fast moving MBPs 
could play a crucial role in AC coronal heating models. In fact, 
fast moving MBPs have a large impact on the amount of energy 
available for AC heating processe s . Thi s was outlined in the the- 
oretical work of Choud huri et al.l d!993l) ; here the authors indi- 
viduated in 'fast motions' (=; 3 kms -1 ) of magnetic footpoints a 
potential source of energy to support the coronal heating. 
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Appendix A: Correction factors for the lifetime 
distribution 

Our true MBPs liftetime distribution is assumed to be N,(t) = 
No exp(-/i ■ t). This distribution is measured with the outlined 
algorithm that has a certain detection probability p to detect a 
feature in a single image. The elements of the principal diagonal 
of table IA.ll give us the fraction of the measured "true" distribu- 
tion in the right measurement bins, given as: 



No • exp(-/t ■ t) ■ p 



(A.l) 



We are now deriving the higher order correction factors for the 
measured distribution. In Table I A. 1 1 we see that there are a lot 
more elements. The first row (true) indicates the number of 
MBPs at a certain time instance t. The first column (measured) 
indicates the measured number of MBPs for a certain time in- 
stance t. Every element of a crossing of a row (refering to a cer- 
tain measured lifetime t m ) with a column (refering to a certain 
true lifetime t t ) gives us the fraction of the number of MBPs with 
a measured lifetime t m but having an actual lifetime of t t (see also 
Table IaTTT i. We see that not only the fraction of features which 
are measured with the right lifetimes contribute to the number 
of MBPs measured at a certain time but also the other elements 
of the row (which live longer as they are actually measured). We 
can conclude that the sum of the row elements contribute to the 
measured number at a certain time. To consider this we introduce 
a correction factor Kl(t): 

N m = No ■ exp(-A ■ t) ■ p-z +l + Kl(f), (A.2) 
where Kl is: 

K1(0) = p • (1 - p) ■ N,(At) +P-0--P)- N,(2 ■ At) 
+p-(l-p)-N,(3-At) + ... 

= P'(l-p)' I N -exp(-A-t)dt 
J At 

Kl(l) = p 2 ■ (1 - p) ■ N(2 ■ At) +p 2 .(l~p). N,(3 ■ At) 
+p 2 -(l-p)-N,(4-At) + ... 

/-*oc 

= p 2 -(l-p)- I N • exp(-/l • t)dt 

JAt-2 

Kl(t) = p-s +l -(l-p)- I N,(0) ■ exp(-A ■ t)dt 

Jt+At 



= No ■ exp(-/l ■ (t + At)) ■ p~' 



A 



We see that this Eq. is very similar to Eq. I A. 1 1 Thus we can 
rewrite Eq. IA.2l as 



N m = No-exp(-A-t)- p~- +1 + Kl(t) 
N m = N ■ exp(-A ■ t) ■ p^ +l • (1 + Kl), 

where 

l-p\ 



Kl 



■ exp(-zlAf) 



(A.3) 
(A.4) 

(A.5) 



is always positive. This value describes the overestimation of 
the measurement bins of a distribution by erroneously measur- 
ing long living features in the short living bins (by loosing track 
of a feature). Interestingly this factor is time-independent and 
a constant for certain parameters (detection probability p, time 
cadence Af and the physical decaying parameter A). 



Now we can derive a second order correction parameter. 
Loosing track of the MBPs does not only lead to erroneously 
measured MBP lifetimes but also creates a second order lifetime 
distribution. This can be explained that by loosing track of the 
MBP, the MBP itself will not disappear. Therefore, it can happen 
that we will measure the same MBP again in the future. We will 
now try to estimate a parameter describing this circumstance. 

If we move from the principle diagonal in Table [ATI to the 
next diagonal, this diagonal would describe the probability to 
have an MBP in the data set that will live for one time step more 
than originally measured. The sum over these elements will give 
us all elements that live one more time step. The next diagonal 
describes the MBPs that live two more time steps and so on. 
Therefore, the diagonal elements describe the truncated part of 
the original distribution which forms itself a new distribution. 
We will derive this distribution now: 



No, 



p ■ (1 - p) ■ No ■ exp(-A ■ At) + 

p 2 ■ (1 - p) ■ No ■ exp(-A ■ 2 ■ At) + . . 

/-»co 

I p^(l - p) ■ No ■ exp(-A- t)dt 

J At 

p ■ (1 - p) ■ No ■ exp{-A ■ 2 ■ At) + 
p 2 ■ (1 - p) ■ No ■ exp(-A ■ 3 ■ At) + . . 

pir l (i - p ) . n q . exp(-A ■ t)dt 

2-At 



f 



> (1 - p) ■ No ■ exp(-/l • r)dr 



f 

= p & '-(l-p)-No I p A ' exp(-/l ■ r)dr 

Jt+At 

pi.(l-p). No ■ exp((-A ■ • (t + At)) 



= K2 ■ No ■ exp(-A ■ t), 
where K2 is given by: 



K2 



1 



In p 

~aT ' ' 



exp(-/lAf + In p). 



(A.6) 
(A.7) 

(A.8) 

(A.9) 
(A. 10) 

(A.ll) 
(A. 12) 
(A.13) 

(A. 14) 
(A.15) 

(A. 16) 



This factor K2 describes the creation of the new distribution 
(the truncated one), which also contributes to the measured dis- 
tribution. This truncated distribution gives rise to a truncated- 
truncated distribution and so forth. Finally, we arrive at the fol- 
lowing measured distribution: 



NJt) = No ■ exp(-A ■ t) ■ p~> +l ■ (1 + Kl) 
•(1 + K2 + K2 2 + K2 3 + . . .) 



y KT = ^— 

,4i . I - K2 



=0; K2<l 



Np-expj-A-t)- pi: +l -(l+Kl) 
l -K2 



(A. 17) 
(A. 18) 

(A. 19) 
(A.20) 



As a last point, we want to give the equation for the cumulated 
number-lifetime measurement distribution of the MBPs: 



Nm cum (t) 



exp(-A-t)- p~- +1 -(1+^1) 
1 -K2 



dt 



1 +KI 



1-K2 \ _ !iL£ 

A At 



(A.21) 
(A.22) 
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Table A.l. shows the probability relations (fractional relation) between the true distribution (N t ; indicated by the first row) and the 
actually measured distribution (N m ; indicated by the first column) 2 . 



measured 



true N(t) 


N,(0) N,(Af) 


N,(2 ■ At) 


N,(3 -At) 


N m (0) 


P p-(l-p) 


p-(l-p) 


P-H-P) ■■■ 


N,„(At) 


P 1 


p'-a-p) 


pM1-p) ... 


N m (2 • At) 






p s -(1-p) ••• 


N,„(3 • At) 









•M,-exp[[-,l,, + ^ ln(p)| /) (A.23) 

- || .1 ( ^ liK/O] /]• (A.24) 

where K cum is the cumulated correction factor: 
_1+*1._P_ 

Kcum ~ l _ K2 l±£ - ( A -^) 

A At 

If we compare this with the true cumulated distribution N,(t) = 
■ exp(-At), where No/ A is an interesting parameter (because it 
describes the total amount of MBPs in the data set), we see that 
we have to divide our measured N mcwn (0) by the factor K cum and 
by A: 

N t (0) = iV " w( ' Q) . (A.26) 



2 As an example, let us think of the real number of MBPs at time 
instance 2At (N,(2At)), this number is measured (contributes) by a fac- 
tor of p 3 reduced for the measured number N m (2At). But not only the 
true number of MBPs at a certain lifetime contribute to the measured 
value of this lifetime. If we consider the measured number of MBPs at 
time At, we see that the true number of MBPs at time 3 • At contributes 
with a fraction of p 2 {\ - p) to the measured number of MBPs at a life- 
time of At. To get a correct distribution all of these relations have to be 
considered (see text). 



